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Abstract: Production of heat shock proteins are induced when a living cell is ex- 
posed to a rise in temperature. The heat shock response of protein DnaK synthesis 
in E.coli for temperature shifts T — > T + AT and T — > T — AT is measured as 
function of the initial temperature T. We observe a reversed heat shock at low 
T. The magnitude of the shock increases when one increase the distance to the 
temperature To ~ 23°, thereby mimicking the non monotous stability of proteins 
at low temperature. This suggest that stability related to hot as well as cold 
unfolding of proteins is directly implemented in the biological control of protein 
folding. 

PACS numbers: 5.65. +b, 82.60-s, 87., 87.10.+e, 87.14.Ee, 87.15.-v, 87.17.-d 



Chaperones direct protein folding in the living cell by binding to unfolded 
or misfolded proteins. The expression level of many of these catalysts of pro- 
tein folding change in response to environmental changes. In particular, 
when a living cell is exposed to a temperature shock the production of these 
proteins are transiently increased. The response is seen for all organisms 
[[I]], and in fact involves related proteins across all biological kingdoms. The 
heat shock response (HS) in E.coli involves about 40 genes that are widely 
dispersed over its genome |2|, [|. For E.coli the response is activated through 
the a 32 protein. The a 32 binds to RNA polymerase (RNAp) where it dis- 
places the a 70 subunit and thereby changes RNAp's affinity to a number of 
promoters in the E.coli genome. This induce production of the heat shock 



proteins. Thus if the gene for a 32 is removed from the E.coli genome the HS 
is suppressed 0, U and also the cell cannot grow above 20°C. 

The HS is fast. In some cases it can be detected by a changed synthesis 
rate of e.g. the chaperone protein DnaK already about a minute after the 
temperature shift. Given that the DnaK protein in itself takes about 45 
seconds to synthesize the observed fast change in DnaK production must be 
very close to the physical mechanism that trigger the response. In fact we will 
argue for a mechanism which does not demand additional synthesis of a 32 in 
spite of the fact that DnaK is only expressed from a a 32 promoter, and thus 
postulate that changed synthesis of a 32 only plays a role in the later stages 
of the HS. To quantify the physical mechanism we measure the dependence 
of HS with initial temperature and find that the magnitude of the shock is 
inversely proportional to the folding stability of a typical globular protein. 



The present work measure the expression of protein DnaK. Steady state 
levels at various temperature and growth conditions can be found in || §]. 
The steady state number of DnaK in an E. coli cell varies from approx. 4000 at 
T = 13.5 to approx. 6000 at 37°C, thereby remaining roughly proportional 
to number of ribosomes. DnaK is a chaperone, and have a high affinity 
for hydrophobic residues J7|, as these signal a possible misfold (for folded 
proteins the hydrophobic residues are typically in the interior), a 32 control 
the expression of DnaK, the a 32 must bind to the RNAp before this bind to 
the promoter for DnaK. One expect at most a few hundred a 32 in the cell, 
a number which is dynamical adjustable because the in vivo half life of a 32 
is short (in steady state it is 0.7 minutes at 42°C and 15 minutes at 22°C 
[0, |T0[]). The life time a 32 is known to increase transiently under the HS. 



The measurement was on E.coli K12 strain grown on A+B medium [jT2 
with a 3 H labelled amino acid added. After the temperature shift we ex- 
tracted 0.6 ml samples of the culture at subsequent times. Each sample was 
exposed to radioactive labelled methionine for 30 seconds, after which non 
radioactive methionine was added in huge excess (10 5 fold). Methionine is 
an amino acid that the bacteria absorb very rapidly, and then use in protein 
synthesis. Protein DnaK was separated by 2-dim gel electrophoresis as de- 



scribed by O'Farrell QTTJ , and the amount of synthesis during the 30 seconds 
of labeled methionine exposure was counted by radioactive activity and nor- 
malized first with respect to 3 H count and then with respect to total protein 
production. This result in an overall accuracy of about 10%. The result is 
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a count of the differential rate of DnaK production (i.e. the fraction DnaK 
constitute of total protein synthesis relative to the same fraction before the 
temperature shift [|I2[]) as function of time after the temperature shift. For 
the shift T — > T + AT at time t = we thus record: 

Rate of DnaK production at time t 
Rate of DnaK production at time t = 

where the denominator counts steady state production of DnaK at temper- 
ature T. In figure 1 we display 3 examples, all associated to temperature 
changes of absolute magnitude AT = 7°C. When changing T from 30°C to 
37°C one observe that r increases to ~ 6 after a time of 0.07 generation. 
Later the expression rate relaxes to normal level again, reflecting that other 
processes counteracts the initial response. When reversing the jump, we see 
the opposite effect, namely a transient decrease in expression rate. Finally 
we also show a temperature jump at a low temperature, and here we observe 
the opposite effect, namely that a T increase give a decrease in expression 
rate. Here a corresponding T decrease in fact gives an increase in expression 
rate (not shown). Thus the cells response to a positive temperature jump is 
opposite at low temperature T than it is at high T. 

In figure 2 we summarize our findings by plotting for a number of positive 
shocks T —>■ T + 7°C the value of r = R where the deviation from r = 1 
is largest. This value can be well fitted by the following dependence on 
temperature T 

ln{R{T)) = (aAT) (T - T ) (2) 

where R(T = T = 23°) = 1 and aAT = iw glff ) = 0.2 • K~ l (i.e. 
a = 0.03K- 2 ). 

In order to interpret this result we first assume that the production rate 
of DnaK is controlled by two factors, a slowly varying factor C that depends 
on composition of some other molecules in the cell, and an instantaneous 
chemical reaction constant K. Thus at time t after a shift in temperature 
the production of DnaK in the cell is: 

< &^M1 (t,T^T + AT) = C{t,T ^T + AT) ■ K(T + AT) (3) 

where the initial composition of molecules, C(t = 0,T — > T + AT) equals 
their equilibrium number at the temperature we changed from, i.e. = C eq (T). 
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To lowest approximation, where we even ignore feedback from changed 
DnaK in the cell until DnaK production rate have reached its peak value: 



_ K(T + AT) 
R ~ K(T) (4) 



which implies that 



ln{R) = ln(K(T + AT)) — ln(K(T)) = AT (5) 

Using the linear approximation in Fig. 2 

ln(K) « const + -(T-T ) 2 (6) 



Identifying K = exp(—AG/T) the effective free energy associated to the 
reaction is 

AG « G - ^(T-T ) 2 (7) 
Thus AG has a maximum at T = T = 23°. 

To interpret the fact that HS is connected to a AG that have a maximum 
at T = T ~ 23° we note that many proteins exhibit a maximum stability at 



T between 10°C and 30°C H, 0. Thus AG = G(folded) - G(un folded) 
connected to the folded state of a protein is at a minimum at T . The corre- 
sponding maximum of stability is in effect the result of a complicated balance 
between destabilization from entropy of polymer degrees of freedom at high 
T, and destabilization due to decreased entropic contribution to hydrophobic 



stabilization of proteins at low T [H , |T5| ]. One should expect a similar be- 



haviour also for some parts of a protein | lq| , and thus expect a max binding 
for hydrophobic protein-protein associations around To. Quantitatively the 
size of the AG change inferred form the measured value of a = 0.03-fT -2 
correspond to a changed G of about 20 — ■> 30kT (about 15Kcal/mol), for 
a temperature shift of about 40-50°G. This matches the change observed 



for typical single domain proteins Thus the HS is associated to a AG 



change equivalent to the destabilization of a typical protein. 

The above picture still leave us with the puzzle that protein binding and 
folding stability is at a maximum around T , whereas the effective AG we 
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observe have a minimum there. This can only be reconciled if the interaction 
we consider is inhibitory. An inhibitory binding that controls the feedback is 
indeed possible 0. To summarize our understanding we in figure 3 display 
the molecular network that we believe is controlling the transient heat shock 
levels of DnaK in the cell. The key inhibitory control mechanism is the 
association of DnaK to a 32 . DnaK binds to unfolded protein residues 0, and 
the amount of DnaK-a 32 association thereby monitor cellular consequences 
of a shift in temperature. 

Impacts of mutants: We have measured the heat shock in a strain where 
the a 32 gene is located on a high copy number plasmid. In this strain where 
the synthesis rate for a 32 may approach that of DnaK we find a HS that 
was smaller and also remained positive down to temperature jumps from 
T well below T = 23. According to fig. 3 this reflects a situation where 
both a 32 and DnaK are increased. With increased DnaK level, one may have 
a situation where DnaK exceeds the amount of unfolded proteins, and free 
DnaK concentration thus becomes nearly independent of the overall state 
of proteins in the cell. Also the huge increase in a 32 supply may decrease 
the possibility for the sink to act effectively. Thereby other effects as e.g. 
the temperature dependence of the binding a 32 ■ RNAp versus the binding 
a 70 -RNAp (i.e. K32/K70 from figure 3) may govern a response that otherwise 
would be masked by a strongly varying inhibition from a 32 binding to DnaK. 



The reaction network in Figure 3 allow a more careful analysis of the 
production rate of DnaK: 



d[DnaK] 



a 



:-!2l 



oc [RNAp ■ a 32 } « ^__J — ( 8 ) 



dt 1 + g [DnaK] 

where the a 32 changes when bound to DnaK due to degradation by proteases 
(the "sink" in figure 3): 



d[a 32 ] 



32 1 



oc Supply - [DnaK -a 32 } « Supply - ° — — r— T (9) 



dt 1 l + (g[DnaK]) 

Here g = exp(—AG/T) is an effective reaction constant. In the approxima- 
tion where we ignore free a 32 , free a 70 and the fraction of DnaK bound by 
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cr 32 then: 

( KroW 70 } \ ( K D32 \ 
9 ~ \K 32 [RNAp]J {l + K DV [U f )J 1 ' 

The first term expresses the cr's competition for RNAp binding whereas the 
second term expresses the DnaK controlled response. [Uf], which denote 
unfolded proteins that are not bound to DnaK, decreases with increasing 
[DnaK]. 

When moving away from T , i.e. lowering g by increasing [Uf], the rate 
for DnaK production increases. For an approximately unchanged "Supply" 
the extremum in production occurs when d[a 32 ]/dt = and has a value that 
approximately is oc l/g. With the assumption that "Supply" does not have 
time to change before extreme response is obtained, we identify R with l/g 
and thereby with the free energy difference AG that controls the HS. The 
early rise in r is reproduced when most cr 32 are bound to RNAp reflected 
in the condition g[DnaK] << 1. This implies a significant increase in a 32 
lifetime under a positive HS, and implies that the early HS is due to a changed 
depletion rate of a 32 . Later the response is modified, partly by a changed 
"Supply" and finally by a changed level of the heat shock induced protease 
HflB that depends on and counteracts the cr 32 level in the cell. 

The largest uncertainty in our analysis is the possibility of a significant 
time variation in "Supply" and HflB level during the HS. As these will govern 
the late stages of the heat shock, in particular including its decline, the 
variation in "AG" for proteins in the cell may easily be underestimated from 
using the peak height variation with T. Adding to the uncertainty in what 
AG precisely represents is also the fact that although we only measure DnaK, 
it can be the complex of the heat shock proteins DnaK, GrpE and DnaJ that 
sense the state of unfolded proteins in the cell, as removal of any of these 
display an increased life time of a 32 |9|, [UJ . Such cooperativity may amplify 
the heat shock. 

For the final interpretation of AG we stress that it effectively counts 
the free energy difference between the complex DnaK-a 32 and that of DnaK 
being free or being bound to unfolded proteins in the cell. Dependent on the 
fraction of DnaK relative to unfolded proteins [U] in the cell, i.e. whether 
Knu\Uf\ is larger or smaller than 1, the HS will depend or not depend on the 
overall folding stability of proteins in the cell. Thus for much more unfolded 
proteins than DnaK in the cell, the measured AG reflect both an increases 
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of the binding to unfolded residues Kou[Uf] as well as a decrease of the 
DnaK-cr 32 binding Kds2 when moving away from To. Our data does not 
discriminate between these processes. This discrimination can however be 
obtained from the data of ref. [16| where it was found that overexpression 



of DnaK through a cr 32 dependent pathway in fact repress HS. As DnaK-cr 32 



binding still play a crucial role in this setup, the vanishing HS of |16| support 
a scenario where too much DnaK imply Kou[Uf] << 1. Then g in eq. 10 
and thereby the a 32 response becomes insensitive to the amount of unfolded 
proteins in the cell. 

We conclude that the HS is induced through the changed folding stability 
of proteins throughout the cell, sensed by a changed need of chaperones. We 
believe this reflect primarily an increased amount of proteins that are on the 
way to become folded (nascent proteins), and not an increased denaturation 
of already folded proteins, because spontaneous denaturation of proteins is 
extremely unlikely at these temperatures. Thus the deduced sensitivity to the 
thermodynamic stability of proteins may primarily reflect a correspondingly 
sensitivity to change in folding times. 

We now discuss related proposals of "cellular thermometers" for the HS. 
McCarthy et al. [I7| proposed that the thermometer was a change in au- 



tophosphorylation of the DnaK protein. This should cause a temperature 
dependent activity of this protein. However their data does not indicate that 
the reversed HS response that we observe at T < 23C could be caused by 
such mechanism. Gross Hj made an extensive network of possible chemical 
feedback mechanisms which connect a rise in cr 32 level with the folding state 
of proteins in the cell. It included HS induced through an increased synthesis 
of cr 32 , an increased release of a 32 from DnaK as well as an increased stabil- 
ity of cr 32 when DnaK gets bound to unfolded protein residues. Our fig. 3 
specifies these possibilities to a minimalistic chemical response including the 
two latter mechanisms combined, and of these only the option of a changed 
stability of a 32 due to a sink controlled by DnaK-cr 32 is able to reproduce also 
the fact that the max HS takes time to develop. In regards to the by Morita 
H| proposed increased synthesis of cr 32 , we note that for high temperatures 



T, the major mechanism that controls cr 32 synthesis in fact is a T dependent 
change in the mRNA structure that leads to an increased translation at in- 
creased T |Tj|. However again our finding of a reversed HS at T < 23C is 
not readily explained by such changes in the stability of mRNA structures 
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below 23C. 



In summa, we observed that positive heat shock is induced when T 
changes away from To ~ 23°. We found that the size of the heat shock 
qualitatively as well as quantitatively follows the thermodynamic stability of 
proteins with temperature. This suggested that stability related to hot as 
well as to cold unfolding of proteins is implemented the in HS. We demon- 
strated that such an implementation was possible in a minimalistic chemical 
network where the control is through an inhibitory binding of the central 
heat shock proteins. We finally saw that the temporal behaviour of the HS is 
reproduced when this inhibitory binding controls the heat shock by exposing 
a 32 to a protease. 
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Figure Captions 

• Fig. I. Heat shock response measures as DnaK rate production change 
as function of time since a temperature shift. The production is nor- 
malized with overall protein production rate, as well as with its initial 
rate. In all cases we use absolute AT = 7°C. We stress that the time 
scale is in units of one bacterial generation measured at the initial tem- 
perature. At T = 37°C the generation time is 50 min, at 30°C it is 75 
min and at 20°C it is 4 hours. 

• Fig. 2. Induction fold R for positive temperature jumps as function of 
initial temperature. The straight line correspond to the fit used in eq. 
2. 
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• Fig. 3. Sufficient molecular network for the early heat shock. All 
dashed lines with arrows in both ends are chemical reactions which 
may reach equilibrium within a few seconds (they represent the home- 
ostatic response). The full directed arrows represent one way reactions, 
with the production of DnaK through the a 32 ■ RNAp complex being 
the central one in this work (this step is catalytic, it involves DNA 
translation etc.) The time and temperature dependence of the early 
HS is reproduced when most DnaK is bound to unfolded proteins, and 
when remaining DnaK binds to a 32 to facilitate a fast depletion of a 32 
through degradation by protease HflB. 
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